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ABSTRACT 


The purpose of this work is to investigate further the 
applicability of the optical model in analysis of very light elastic 
scattering systems. 

The optical model was used to fit the elastic scattering 
problems: neutron-tritium from 6.0 MeV to 23.0 MeV, and the proton- 
helium-3 problem from 5.51 MeV to 156 MeV. 

The analysis of the data yielded fits to the cross section 
and polarization data which were usually as good as or better than 
those obtained by phase shift analysis or the resonating group method. 
Beyond the use of the standard optical model, exchange terms were 
added and used in a reanalysis of the p +3He scattering data. The 
addition of the exchange terms improved the backward angle fits 
markedly at 30 MeV and 156 MeV. 

It has been shown here that the optical model is a valid 
method in the analysis of light nuclei scattering and is in many cases 


easier to use. 
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CHAPTER 1 
INTRODUCTION 


The optical model has generally been thought to apply only to 
nucleon scattering on medium to heavy nuclei. In heavy nuclei, the 
energy levels are more closely spaced and therefore, for nucleons with 
energy = 10 MeV, there are no isolated sharp resonances to affect the 
interaction and the process is describable in terms of a complex 
potential. Recently, however, several authors (Sa68, De72) have used 
the optical model to analyze scattering data for very light nuclei with 
a fair amount of success. 


Su, Aue, 


The usual methods of analysis of light nuclei, e.g. She, 
are phase shift analyses and resonating group calculations. Phase shift 
analysis is essentially a phenomenological analysis of data determining 
the phase shifts for a series of partial waves. This method tends to be 
cumbersome since a large number of phase shifts and mixing parameters 
are needed and many extraneous solution sets are found. The optical 
model is much simpler to use, since once the well geometry is determined, 
at most five parameters are searched over and in most cases a single 
solution is found. The resonating group calculations are better in 
principle since they essentially start from the fundamental two-body 
interaction. Generally, however, only one channel is considered in the 
analysis (two channels have recently been used in Mc74 and Ra74) 
ignoring all other channels. Also, when analyses are made at’ higher 
energies, a phenomonological imaginary potential is introduced in order 
to take into account the effects of the other open channels. Thus, the 
advantage of the resonating group calculations being microscopic and 


essentially parameter-free is partially lost. The optical model can be 
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Z 
altered to incorporate most of the effects included in the theoretical 


calculations. In addition, there is another advantage in using the 
optical model. In reactions such as (p,n), (p,2p), and (p,pd), the 
distorted wave Born approximation is used in analysing the data and the 
optical model parameters are needed in order to determine the distorted 
waves (Ma74). 

The present work consists of the optical model analyses of 
almost all of the available data for the elastic scattering systems 


n +34 and p re 


He. The resonances that have been "observed" for these 
systems (1065, To66) are very wide, therefore, it is expected that the 
model is applicable. 

Chapters 2 and 3 are concerned with a brief discussion of 
phase shift analysis and the resonating group method (and other 
microscopic calculations), respectively. The theory of the optical 
model is presented in Chapter 4, followed by the optical model analyses 
of the n 4 and p He elastic scattering problems in Chapter 5. In 
Chapter 6, a discussion of exchange effects is given. These are then 
incorporated within the optical model framework. Chapter 7 entails a 
discussion of the results of the optical model analysis as compared 


with those obtained using the other methods of calculation, followed by 


the concluding comments. 
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CHAPTER 2 
PHASE SHIFT ANALYSIS 


One of the most common methods of data analysis for scatter- 
ing off light nuclei is phase shift analysis. A wave function can be 
expanded into partial waves. A phase shift is essentially the phase 
angle through which the scattered partial wave differs from the 
unscattered one at large distances from the scattering source. Hence, 
if one knows the phase shifts for all of the partial waves, one can 
describe the scattering exactly. 

As a simple example, one can consider a spin-zero scattering 
from a spherical potential V(r). Expanding in a series of Legendre 


polynomials, we have the wave: 


v(r.e) = J ae P, (cosé) (2501) 
R 


and the scattering amplitude: 


Fo) en) FP, (cose) (2.2) 
£ 


Let 
ey ee all fl (2.3) 


now yy(r) is a regular solution of the radial equation: 


2 
d 2 (R41 
(k* - U(r) - = 0 (2.4) 
{> + r , } yp(r) 


and it has the asymptotic form: 
. ] 
‘igeoras sin(kr - 5 21+ 5.) (2.5) 
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4 
~ where So is the phase shift and ao is a normalization Pe eeani: We 
can now write the asymptotic form of w as a series of Legendre poly- 
nomials: 

ik-r ikr 


p(rs8) =e MSS FG) 


ik 
J (200) 5, (kr) 4 is fp, (cose) : (2.6) 


Now , substituting in the asymptotic form of J, (kr) and separating 


the incoming and outgoing waves, one gets: 
rv(r,6) x yf (1) oat apake 
Q 
+ (Bl FY elk" bp, (cos) | (2.7) 


The asymptotic form of Yo equals the bracketed quantity on the right 


hand side of Equation (2.7), therefore, one has: 


ié 
eet, Oeele ee 
gbaber re: (2.8) 
and 
is 
f= aut! sinters. (2.9) 


Substituting these into Equation (2.2) one obtains f(@) as a function 
of 60: 

el=E J (2u+1)e sing, P, (cose). (2.10) 
The differential cross section is equal to: 


Se = 1F(0)I° 
(cose )P, . (cose) 
(2ett) 
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- For an alternate development, refer to Appendix I. 

Phase shift analyses have been performed for p + He elastic 
scattering for incident proton energies from 1 MeV to 50 MeV (To65, 
Ha70, Mc70, Ba71). These included the experimental differential cross 
section data as well as their experimental and extrapolated polariza- 
tion values in the analyses. Table 1 gives the phase shift values 
obtained by Harbison et al (Ha70) for the 30 MeV data and the resul- 
tant curves are shown in Figure 1. The type 1 solutions are those 
which use as initial parameters, phase shifts calculated by the 
resonating group method. This solution set yields s-wave phase shifts 
which are consistent with those found at very low energies. 

An analysis of the n+ 3H problem was performed by Tombrello 
(To66) for 6 MeV neutrons. Parameters from the pt She analysis were 
used to begin the search. The phase shifts obtained are given in 
Table 1 and the fit to the differential cross section is shown in 
Figure 2. 

Trombrello also performed single-level analyses for p+ “He 
(To65) and n+3H (To66) problems using a hard core radius of 4.0 fm. 
The resonance parameters are given in Table 2. 

One problem with the phase shift analysis method is that 
there are upwards of eleven parameters to be searched over, including 
phase shifts and mixing parameters, even restricting the analysis to 
g<2. (The channel-spin-mixing parameters e« account for coupling 
between singlet and triplet states having the same L and J values 

- pt 


such that E79. 37 J” = 17, 2°, 3°. The tensor couplings of states 


with the same J", but different L values are described by the 
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TABLE 1 6 
TYPE I PHASE SHIFTS FOR E, = 30.5 MeV AND FOR E, » 6 MeV. 
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p +3He SOLUTIONS 


N +H 
63, 1(A) 1(B) 1(C) 1(D) 1(E) 
Rc URIOLES® | cISe Oe =117.5° —-189.0° S116.097) +00. 7e 
Bia =12310° 423.00 2126550 4 -W16007 =1133.0°. Leja 29 
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TABLE 2 
RESONANCE PARAMETERS 


SYSTEM CHANNEL ANGULAR CENTER OF MASS RESONANCE 


SPIN MOMENTUM RESONANCE WIDTH 

. ENERGY 42 

S ‘) (MeV) (MeV ) 

p +7He 1 a 47.4 5.5 
1 1" 6015 a5 

1 07 7.94 ane 

0 is 9.79 8.8 

aceon oe 3.4 5.5 
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Figure 2 Phase shift fit to 6 MeV n +4 differential cross 
section. Reference (T066). 
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~ coupling parameters Too» connecting 75 and cba and 113° connecting 
°P, and ae) 

A major problem is that often many sets of solutions are 
found which yield equally acceptable fits, for example, there were 
five type 1 solutions found for the 30 MeV and 49.5 MeV (Ha70) cases. 
Until further data is available for analysis the extraneous solutions 
cannot be eliminated. Having several solutions leads to difficulty 
in extracting reliable physical significance from the resultant phase 


shifts. 
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CHAPTER 3 
RESONATING GROUP METHOD 


The method of resonating group structure has been used to 
calculate cross sections (polarizations are not done) for scattering 
systems involving light nuclei. It essentially consists of building 
a completely antisymmetrized wave function for the whole system, 
target plus projectile, from channel functions which correspond to 
the various possible ways of grouping the nucleons, for example, oa, 
She, d, etc. (Ta63). This method investigates the close interactions 
within the nucleus between these groups as opposed to looking at 
average potentials. 

The wave function py for a total system of m protons and n 
neutrons (m+n=N) is written as a sum of terms, where each term 
represents the N nucleons grouped into a particular configuration. 
Each term is a product of the wave function ¢ describing the motion 
of the particles within each group, multiplied by a function au 
depending on the positions and spin variables Mm. (total internal 
angular momentum) of the different groups. The properly anti- 
symmetrized wave function for the whole system, assuming the groups 
are spinless, has the form: 


wal ] 
w(1,2,...N)=a, } (+1) d, F °Lm, *I1 4m 
m 


perm Il * 


SRE Nig Ege ek Peg We len 
perm me tag 


where » (+1) is a summation over all m! permutations of the pro- 
pérm 


tons and n! permutations of the neutrons with a change in sign for 
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odd permutations; the a; are factors simplifying the normalization; 
the oF om. are the different states of the same group and are ortho- 
normal , the m, singles out the wave function for group j which 
represents the z component of its angular momentum as mss the d; 
gives a partial wave function wy corresponding to definite values of 
the angular momentum of the whole system and its projection along the 
z-axis. 

The main problem is the calculation of the paige They are 
determined uniquely by demanding that they yield the best wave func- 
tion w, with respect to the variational principle: 


6E = 0 where Bots LUxH Vdt : (3.2) 
fuxwdt 


This yields a set of simultaneous integro-differential equations in 
the ee which can be transformed into Fredholm integral equations 
through the use of a generalized Green's function. The condition 
that the Fredholm determinant for this set of equations vanishes 
determines the phase shifts, as well as energy levels and transmuta- 
tion probabilities. 

The integro-differential equation is obtained from the 


Rayleigh-Ritz variational principle in the following way. We have: 


sate 5 [ue (H- Ean (3.3) 
now the wave function is of the form: 
v = o(r)F(R) (3.4) 


where is in internal (relative) coordinates and F is in external 


(center of mass) coordinates. 
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Before carrying the treatment any further, let us for sim- 
plicity, take the p + He scattering as an example. To include the 
exchange of the incident proton and one of the She protons, the wave 


function must include terms of the form: 
6(1,2,3)F(0,°He) - 6(0,2,3)F(1, He) (3.5) 


where the Sie contains nucleons 1, 2, 3 and the incident particle is 
0; 0, 1, and 2 are protons and 3 is a neutron. The wave function is 
itself antisymmetrized with respect to the two protons. 


Allowing for full exchange the wave function becomes: 
v = $(1,2,3)F(0, He) - 9(0,2,3)F (1, He) +9(0,1,3)F(2, He) . (3.6) 


To indicate how the exact calculation is done, a simpler derivation 
of the integro-differential equation will follow using the simpler 


wave function of Equation (3.5). For convenience let: 


1 = (12,3) ge =a(0j2}8) 
(37) 
ia = F(0,>He) F, = (1, He) 
9) ] 
therefore, 
eS Fy = doy . (3.8) 
We have the Hamiltonian: 
= li eee 3.9 
Hs he tah ; (9) 


where ho is the internal hamiltonian for the projectile which we 
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may leave out since it is a single particle; h, is the internal 


hamiltonian for the target: 


3 
j ES V6 
icf te 


such that h,o = ey(we shall however set «=0 for convenience), 

H' is the interaction hamiltonian between the target and projectile 
where 
SA ol 

to is the kinetic energy of the relative motion between the projectile 


and the target: 


2 ne 
“Om ’ 


t. is the kinetic energy of particle i: 


cae? 
em 5 


and ee is the interaction between particles i and j. 


Therefore, the Hamiltonian may be written: 


: eee 
H = om V = yy (3.10) 
where 
= 5 A 3h 
velv, (3.11) 


j 
Substituting the above equations into Equation (3.3), one obtains 


the equation: 


* * 
6E 6 (OF o- oF) (H- E) (oF, - oF) dP 553 
r 6] (Foe) - FTPQ)H (gf - O4FQ) at - E8| (Fao - F700) (oF 1- oy Fo) oss 


(3.12) 
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Since an equation in the external coordinates is required, the varia- 


° ° * 
tion is performed over FO therefore, 


2 
_ fh gaz * omen e® 
SE = Fy | teFev obo oFe|6TVO Fy) 


Simos 

F o> oro or" 6 Fy] +| 19F 
* el ne 3 

- E| (OFSF, - FS | 046 Fy] (3.48) 


where 
V. = 1o*Vo_ d 
Di 2) one ce ie3 


is the average potential over the internal structure. Now, for 


sE =0, 
f° RA eer aie ode, WiC Ee Uae lec V sk: dra cl 
om | orrvo 11'6 4 %% qons3 OuD.o8t. OFT %% 1~ 133 
* xl a* a 
Saks aa oF4 [oto Fydry al mht ss (3.14) 
Therefore, 
else? Mee meee ey BENG F.dt (3.15) 
2m Dir rf) 12m 0 1 ; : 
r' 
Now 
(fe yy Ey6u= Kiron) = kernel (3-16) 
Ooms - . = Ger = kerne 5 


where r is for the unexchanged case and r' is for the exchanged case. 


Therefore, 


citngting ee [k(r rE (r')dr! (3.17) 
nD r ror! )F(r 


~ 


which is the integro-differential equation we are looking for. This 
equation is not usually used in this form, the functions in Equation 


(3.17) are expanded in terms of partial waves. We write: 
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2 2 2 2 2 
-f 2h. ~ e Ac a , ATR(L+1 
om Y= Oi earner a (3.18) 
r 2mr 
f(r) 
Pir) => ) P. (cos 6) (3.19) 
hy q r R 
and 
Ka bner) 
Kors Tee m* Mot 
sali ight" til Yo (2)¥ 5 (2 ) (3%.20) 
Therefore, we can write Equation (3.17) as: 
NEE! 0% She Hawa) 
) ~~ —3+ Bat Vee E]f, (r)P, (cos 8) 
Q or 2mr 
Kr we) 
EP ee Clg Qi? m* Moi ; 
a‘ ! | I ar YF (2)N(2')P, (cos!) (3.21) 


Equate terms in & and drop the 1/r, then left multiply both sides by 


* 
y°"(2) and integrate over 2. Recall that: 


Q 

P,(cosé) = [ga ¥9(2) (3.22) 
also 

m*,m' 7 

\v" Vy 25 5p 1G (3.23) 
and 

ym (2) = (-1)"5"\(2) (3,24) 
Therefore: 

De 2 

ee ae Ett (r) = [k (rye!) (r!)dr' (3.25) 
ar ee are D Q Lae BAS, 


Resonating group calculations in the one-channel approxima- 
tion have been done for p+ "He and n+ 34 elastic scattering systems 
by Reichstein, Thompson and Tang (Re71). They used two different 
internal wave functions for the She, a single gaussian and a two 


gaussian: 
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1tha® 2 
$= Expl- see (re-R)“] (3.26) 
= 
and 
1 3 2 1 3 2 
a koa ica 1ectbe, Gea Puts 2-5 “ideo lee bs (3.27) 
i= i= 


where the a's are variables, the in is the position vector of the jun 


nucleon from the projectile and R is the position vector of the center 
of mass of the three nucleon cluster. They used a nucleon-nucleon 


potential of the form: 


ide Oo 1 Oo Ties lye. r 
MG ge Qt HPs WE TOURS VST (put {2 u)P; 5] 
ef 
+ ary (1 +125) + Ty) ‘tn (Shee) 


where PS 5 is the spin-exchange operator, PY. is the space-exchange 


ij 
operator and Tae is the z-component of the isospin operator for the 


jen particle. The quantities Vy and ve are the s-state triplet and 


Singlet potentials given by 


V.=-Vie (3.29) 


Vo=- Vo e . (3.30) 


where the constants were found to be (Re71): 


2 


V5 = 66.92 Mey Ky 0.415 fm 


ot 


0.292 fm-2 


Vee 29.05 MeV Ke 


and finally, u is an adjustable parameter for best fit. 
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They found that the two gaussian wave function gave slightly 


better fits, in general. They also noted that the effective triplet- 
state p + He potentials in the odd-& states are quite different from | 
those in the even-2 states, see Figure 3. Some of the fits to the 
experimental p+ *He data they found are shown in Figure 4. The n+ 34 
data were analysed in a later paper by Le Mere et al (Le75) and some 
of their curves are shown in Figure 5. The phase shifts calculated 
using this method are useful as initial parameters for phase shift 
analysis (see Chapter 1). The phase shifts obtained by Reichstein 
et al for 30.5 MeV protons from Site are compared in Table 3 with a 
set obtained by Harbison et al (Ha70). 

The resonating group method has several advantages in that 
it employs a nucleon-nucleon potential, as opposed to a smooth 
projectile-target average potential. This potential includes exchange 
terms which have been shown to be important (Re71) and these include 
the Pauli principle accurately since the wave function is completely 
antisymmetric. However, due to computational difficulties, approxi- 
mations are made, such as considering only a single channel, and a 
relatively simple nucleon-nucleon potential, which cause problems with 
accuracy, particularly at higher energies. Also, since the nucleon- 
nucleon potential is purely central, phase shift splitting, channel- 


spin mixing and polarization cannot be obtained. 


OTHER MICROSCOPIC CALCULATIONS 


A further refinement to these calculations consists of 
including all two-body channels, open and closed, in the calculation. 


Hackenbroich and Heiss have performed such calculations for the 
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Figure 3 Effective triplet-state p He potentials 
calculated at 6 MeV with the two guassian 
wave function (solid lines) and the one 
guassian wave function (dashed lines) for 
1 = 0 to 3. Reference (Re71). 
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Figure 5 Comparison of calculated cross section curves for n +9 


H 
scattering with the experimental data points at 19.5 MeV 
21 MeV and 23 MeV. Reference (Le75). 
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TABLE 3 


Phase shifts obtained for 30.5 MeV p +3 


He 
elastic scattering by phase shift analysis 
(Ha70) and resonating group method (Re71). 
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23 
n+3H (Ha71) and p+ He (He72) elastic scattering systems. 


For each channel included in the calculation, a coupled 
integro-differential equation comes out of the Rayleigh-Schrodinger 
variational principle. The resulting set of coupled equations is, 
in general, unsolvable. Hackenbroich and Heiss devised an approxi- 
mation which enables them to solve the set of equations (He70). 
They refer to this methad as microscopic calculations. 

For n open channels, one gets n linearly independent solu- 


tions xy which satisfy the following boundary conditions: 


L; +] 
R. elas )~ Rs for R. > 0 (3.31) 
i Sheth pu(RE) ck ag (R.) for open channels 
x2(R, > ij L. i i L. 1 R, +o (3,32) 
0 for closed channels 


They chose a trial wave function of the form: 


$..f, tad Ja + i al A for open channels 
Pe ies) ot (3.33) 


qq e = ach itp for closed channels 
we] Ju 


where life is a regular coulomb function, gp is an irregular coulomb 
function, d, is a regularized coulomb function and 
the ¥,, are trial functions of the form: 
ees 2 
i = R, exp(-8 Ry) ; (3534) 
The total variational function for boundary condition j, for n-open 


channels and M trial wave functions, is: 
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24 
where 


ss fas Ath. f (3.36) 


}; 
Lj 


G. i ae (3037) 


A is the antisymmetry operator and A. represents the angular momentum 
and internal motion of the clusters. The coefficients a are related 
to the transition matrix elements between the channels i and j (for 
example, for elastic scattering, we can express the ay's in terms of 
phase shifts.) This function is then used in the Kohn-Hulthéns varia- 


tional equation: 
6<w|H- E|y> = 0. (3.38) 


In particular, for the n+ 3H and p+ eHe, for a given vie 


their ansatz is as follows: 


L L 
v =AL) ote, ns ) bene Le (3.39) 
r u 
where 4, 6 are the usual surface functions (fixed during the scatter- 
ing process) and X), and hy are the radial parts of the relative 


motion functions. For example, for aie 


J.T . 3 A 
o2e™ = J C(LS.dsmme)o(FHe)e, 4, ¥, 22) (3.40) 
my > eee aaa | 
and for 4u. 
JeT _ ¥ C(L.S.J,m, me) Guyy (a) (3.41) 
Ser ) Pees. 1 
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function with channel spin S, and p is the translationally invariant 


3 3 


spatial part of the ~He and ~H wave functions. Now p was chosen to. 


be: 
5 ae 
Oi = AC OXD (20. Fy Or (3.42) 
pete Re oe, 
where for She: 
a, = 0.2042 by = 0.04 fm” 2 
a, = 0.6938 b, = 0.11 fm7¢ 
2 = 0. 2 
and for 3. 
a, = 0.02 b, = 0.04 fm“ 
an = 0207168 Bee Olli tem 6 
2 ° 2 e 


Heiss and Hackenbroich made a special choice for the central 
part of the nucleon-nucleon potential, i.e. a local soft core poten- 
tial composed of a superposition of three gaussians, one repulsive 
and two attractive (Ei66), 

: 2 
V(r) = aX a, exp(-b r ) (3.43) 


where a sample set of parameters is: 


a, = 880 MeV b, = 5.40 fm” © 
“ -2 

ay = -70 MeV bo = 0.64 fm 
‘i f -2 

a3 = -21 MeV b3 = 0.48 fm 


With these very light nuclei they also use the matrix elements of the 


simplified potentials with the simplified wave functions (pure s-waves 


* 
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for A < 5). 
The repulsive soft core contained in the central potential 
necessitated the inclusion of short-range correlations in their wave 


functions. They chose these to be of the Jastrow type: 


Bem tr -ri i (3.44) 
mon mon 


where 
f(r) = 1 - a exp(-br) (3.45) 


and _2 
a=70,6,  br=-3.0° fm 


This was approximated in order to shift the short range correlations 


V2 matrix elements 


to the Hamiltonian FHF > Hor? since they obtained 2 
otherwise. 

The fits for the p+ He problem are shown in Figure 6 for 
8.82 MeV and 10.77 MeV. Both the cross section and polarization fits 
are excellent, in general, with the cross section fits slightly too 
small in the backward angles. The differential cross section fits 
for n+3H at 6.0 MeV and 10.5 MeV are shown in Figure 7 and reproduce 
the data quite well from above 80° to 180°, being too small in the 
forward direction. 

The microscopic calculations differ from the resonating 


group calculations in that they assume a form for the wave function 


and all two-body channels are included specifically. 
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Differential cross section (a) and polarization (b) 


fits for p +-He at 5.51 MeV, 6.82 MeV, 8.82 MeV and 


10.77 MeV using microscopic calculation. Reference 


(He72). 
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180 


H at 6 MeV and 14 MeV 


using microscopic calculation. Refence (Ha71). 


CHAPTER 4 
THEORETICAL DEVELOPMENT OF THE OPTICAL MODEL 


The semiclassical optical model was first proposed by Serber 
in 1947 (Se47). This model is referred to as the optical model 
because of the similarity between the scattering and absorption of 
particles by a nucleus and the scattering and absorption of light 
by a cloudy crystal ball. 


Consider a wave of wave number 
2 deme 
k = a (4.1) 


where E is the incident energy and m is the mass of the projectile, 
incident on a nucleus represented by a complex square well of depth 
V and radius R. The particle moves more rapidly inside the nucleus 


due to the attractive nuclear field; it has a new wave number: 


[Raa eee V4 | 
ae (4,2) 


Since the potential V is complex, so is the wave number. 


The refractive index of nuclear matter is given by: 


ie a (4.3) 
iy k Ok : 
At high energies, n x 1 and we can write 

n® = 2(n-1) +1 (4.4) 


since..kR, >>. 1, ky << k and D << k, if we let 
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V=U+ iW (4.5) 
then 
U = vhk, and We 5 VAD (4.6) 
2h 
Vir Th (4.7) 


where v is the velocity of the incident particle. 


Therefore, the nucleon wave function inside the nucleus is: 


a einkz (4.8) 
i(ktky)z -5 
=e e : (4.9) 
From Equation (4.6), one can see that the nucleon wave is attenuated 
by: 
lv.|“ =e : (4.10) 


Using this model, absorption and elastic scattering cross sections 


for neutrons scattered by a nucleus of radius R, are given by (Ho71): 


oy = mL! + [1 - see exp(-2DR i} (4.11) 
2D°R 
R 
og =2n | |1 = exp{(-D+2ik, )S1|? dp (4.12) 
0 
where 
52 ebperet ofl (4.13) 


Scattering by a complex potential can be treated quantum 
mechanically as well by allowing the potential in the Schrodinger 


equation to be complex: 


Vir) = U(r) HINGE). (4.14) 


3] 
A complex potential implies that the wave functions u(r), where 


u(r) 
We ) ine e ey LeOSi0)" 5 (4.15) 
L 


scattering matrix elements S} > phase shifts 6) > and scattering ampli- 
tudes f(0) are all complex as well. 


The Schrodinger equation is of the form 
vp +O (E-U-iWy=0 . (4.16) 
h 


(Note, E is the total real energy of the system, see Equation (4.1).) 
Multiply Equation (4.16) by w* and subtract w times the complex con- 


jugate of (4.16) to obtain: 


vtvey - wey = 2 wr. (4.17) 


Now, recalling the quantum mechanical expression for particle flux j 
and density p, one sees that Equation (4.17) is equivalent to the 


classical equation: 


a div j = Dvo (4.18) 


where v is the velocity of the projectile inside the nucleus. This 
is the well known continuity equation with a sink term on the right 
hand side. Hence, for W(r)<0O, the imaginary part of the potential 


effectively absorbs flux from the beam. 


FESHBACH THEORY 
Consider the Schrodinger equation 
Ey = Hy (4.19) 


where H= He V, and Ho (r) is the unperturbed nuclear Hamiltonian, 
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32 
and V is the projectile-target interaction. 


The total wave function can be expanded in terms of the 


scattering states w of the projectile and eigenstates of Hy? 
wlEsr) = J 6, (Ev, (r) 
- x 


where € are the internal coordinates of the target, and r is the 
relative coordinate between the nucleon and the center of mass of 
the target. Next introduce the projection operators P and Q, 
where 


peer k = -Q (4.20) 


and P projects out the ground state part of the target wave function. 


The Schrodinger equation reduces to the coupled equations: 


(E - Hog)Qv = Hop Py (4.22) 
where 
Hpq = DOO s OLeo. (4.23) 


Eliminating Qy and applying Equation (4.20), yields an equation for 
elastic scattering of the projectile: 


=] 
{E~ H-<b glV[9,> - <¢g|VQ(E- Hog)” QV 90>} Vy = 0 (4.24) 


where it is assumed that the energy of the nuclear ground state is 


ZEO. 


Therefore, the generalized optical potential is given by 


-1 
m3 <ooVIo9> + <$, [VQ(E- Hog) QV 95> (4.25) 
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3S 
where the second term varies rapidly with energy when E is near an 


eigenvalue of Hog giving rise to a resonance in the cross section. 
The propagator (E - fe can be expanded in eigenstates 
of Hag? hence, for scattering near an isolated resonance (E= E.), 


the optical potential is of the form: 


<b WAIer<VelW 140? <rglVAlY,><¥, |QV19,> 


Vopt 4 SO Aa p ah E - E. E - E. 
(4.26) 
<b, [VQ]b,><¥, | QV] o,> 
= <$5/V']¢,2> + = = (4.27) 
S 


where V' is chosen such that <oo1V'1$5> is equal to the first two 
terms in the right hand side of Equation (4.26), and its correspond- 
ing cross section varies smoothly with energy and includes the effects 
of distant resonances. 

Let Yo be the solution to the Schrodinger equation with V' 
as the potential, and let i be the plane wave solution. From the 


Gellmann and Goldberger relation, the transition matrix is: 


, W 

Ti cuba le Ohms aes CVn a ae 
where 

<p |VQ|$.><$_ |QV|o,> 
= on eee dae (4.29) 
S 

and 

@ ture)! s, (4,30) 


Now T can be expressed in Breit-Wigner form: 
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+ Woh ol VOI o,><4 | W1¢ bor 


T=T (4.31) 
O Pats 
Bots Mats! fg 
where 
P 
Ae * <o51W [96> ETH WT “bol YO o,2 a ay 
P, = 2m] Wo >5(E- T- V')<o, | VQ] o.> (4.33) 


P is the principle value and 14 is the scattering amplitude due to V' 
alone and is given by the first term on the right hand side of 
Equation (4.28). 

To derive an optical potential from Equation (4.27), it is 
necessary to average over the energy using an interval A which is 
large compared to the average resonance spacing. The optical model 
potential is defined such that its T-matrix corresponds to the energy 
average of the actual T-matrix. Therefore, average T using a nor- 


malized wave function o(E,E'): 
ae [o(E,E")T(E")ae: (4.34) 


Choose the Lorentzian form for 


1/20 
De) + lee (4.35) 
‘epg (Es) 2(172)- 


where I x 2A/n, and from Cauchy's theorem, we have the relation for 


the weight function: 
fo(E,E )F(E')de* = F(E-3 1) (4.36) 


Therefore, we have: 
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35 
VEE = V+ iW 


(bo VQ19,)(9,] WV 4,) 


V(r) =(o,|V'}o,) + (4.37) 
: : az : E-E.+y iI 
Siz 
where C is the position of the incident particle. 
The imaginary part of the potential is given by: 
2. (bglVOlo,)(o,1QV]o,) 
W= - rh) Fh ee Ree woe 6. . (4.38) 
S Ut ALE or ae) /I 


Now the denominator damps contributions from al] levels such that 


|E-E.| >> 1/2, therefore: 


N 
WF J Wol¥O18, (9.10104) (4,39) 


In practice, however, this expression is difficult to evaluate. 


GREENLEES THEORY 


The interaction between the target and the projectile is 
the sum of all the projectile-nucleon-target-nucleon interactions 
AGS eal?) (See Figure 8). 

Look at the first term in Equation (4.26): 


now, 
A 
Vie) Meer) (4.41) 
i=] peer 
therefore, 


A 
Vopt ~ <b L, vines) o> (4.42) 
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Figure 8 Diagram illustrating vectors used in theoretical 


derivation. 
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37 
or 


Vir) = fog(e')v(Ir= [ar (4.43) 


where mae) is the density distribution of nuclear matter. 

Now, V(r) is the first term in the Feshbach expansion. Greenlees 

et al (Gr68, Gr70a, Gr70b) ignored all the correction terms and used 
Equation (4.43)as their potential. This treatment leads only to the 
real parts of the optical potential. They then used the following 


form for the nucleon-nucleon interaction: 


V(r) 


~ 


ay Cale 


= Vy lr) eV (r)ts tet V r)gy os #V5 (rog-og)(t4 +t) 


+ VO) Vi (edt 415; ; 


wl a 


] 
* a V(r) rye y)* (Ps -py) (a, 404) ] (4.44) 
where 34 is the tensor force operator and 


S54 = 3(0, +) (a5-1) - 05°9; (4.45) 


also: o's are the Pauli spin operators, t's are the isospin opera- 


tors, and P refers to the particle momentum operator. For nuclei 


of spin zero, one can write the potential in the form: 


V(r) = Vor) + Vir) + Vo (r) (4.46) 


~ ~ 


where Vp is purely central, and Vy and Ve contain spin and isospin 


terms. 
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To find expressions for these potentials, breakdown the nuclear 


matter density into its proton and neutron components: 
Pat) = ep(r) + ontr) (4.47) 


Assume spherical distributions and substitute Equations (4.44) and 


(4.47) into Equation (4.43). This gives us: 


Valt) = [og(r' Wp (Ir-r' ar (4.48) 


Vp(r) = 1] top(r')~ or" )IV, (ren! | ar (4.49) 


©0 2n-3 
V(r) u -[f } At én c2(n+1) d (2) 


nay (2neT): rg en-3 Om 
2n-2 'P 

ns aad Pf" )} | 4 s(n)ne*2an |g : (4.50) 
Yr 


Under 50 MeV, it is a good approximation to use only the first term 


for Ve which gives: 
Ver) = og oe ee fr) Ve Cadnian LG (4.51) 
S K 3 ry dr °m S 7 ae : 


Forms for the nuclear distributions and the radial functions 


are needed to perform calculations. These are usually chosen as 


follows: 
alae 
p»(r) = p,(0)/1 + exp( a ) (4.52) 
-ur 
Vir) = Vp e = (4.53) 


where Vp and wp are adjustable. Also assume: 
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op(r) = + 0, (r) (4.54) 
o(r) = * 5 (r) (4.55) 
n A *m ‘ 
(hr eee ae (4.56) 
Z A 7 eR : 
if 
VAr) = ove (r) 
+] for neutrons 
and vas 
-| for protons 


In their analysis of 30 MeV proton scattering experiments by a range 
of nuclei, Greenlees et al (Gr68) used Tm? @m? Vp» Wo» Wy i and a, 
as variable parameters and obtained fits to the data as successful 


as those obtained by other methods. 
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CHAPTER 5 
ANALYSIS OF THE EXPERIMENTAL DATA 


Optical model analyses of elastic scattering by light nuclei 
have been performed in the past with varying degrees of success. 


The potentials used were generally of the form: 


" ; d 
Utr) ¥2 “Vof(rsR,2a,) - iW F(R. sa, ) +4ia:Wp an f(r,R. a, ) 


fid,2 1. “d 
+ no) r or F(r Roa ork (5.4) 
where 
r-Ry -] 
F(r.R, sa) = {1 + exp( i )} (522) 
and 
7 es 
Ry = by A : (5.3) 


In 1964, Kim et al (Ki64) analysed the differential cross section of 
p+ He at 31 MeV without much success. In their potential they chose 
Ro = R. = Reo? and a surface peaked gaussian was used instead of the 
Woods-Saxon derivative. In 1968, n+ “He and p+ “He data were analysed 
by Satchler et al (Sa68) for energies from 0.3 MeV to 18 MeV, and 

fits were obtained which closely reproduced the experimental data. 

The energies considered were below the threshold energy for tHe 
breakup, therefore, there were no imaginary terms in the potential. 
In 1970, Thompson et al (Th70) analysed p + He at energies between 
31 MeV and 55 MeV successfully. Recently, Devries et al (De72) 
analysed p+d using a spin independent optical potential but they 


had to include exchange effects to obtain good fits. They also 
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4] 
attempted to fit the n+ 3H 18 MeV cross section data. 


34 and p + He elastic 


In the present analysis of the n+ 
scattering data, an optical model of the above form is used. The 
Parameters were obtained using the search code SNOOPY (Sc69) and 
were refined with the double precision code TASHTASH (Shpc). The 
final parameters were determined generally by minimum values for 


chi-squared values: 


x2 x2 
va = i + rr (5.4) 
o Pp 


where No and No are the numbers of data points for differential cross 


section and polarization, also: 


OMG.) -oeeyh © 
ex: 7 Esai | (5.5) 


“fe = vo 
7 421 Ad (0; ) 


and O,,(0;), O44, (9, ), and Ao, (85) are the experimental, theoretical 
and experimental error values of the cross section at a center of 
mass angle 0. similarly for Ke. 
The analysis was begun with the n+ 34 problem (Sh72). The 
polarization data at 22.1 MeV was analysed first using initial 
parameters from Devries et al (De72) and estimated values for the 
spin-orbit parameters. The well geometry was determined and then 
fixed for ALL further calculations. Next, the well parameters were 
determined and used as the initial values for the neighbouring 
energies. All available data was analysed. In the p+ He (Po74) 


analysis, the searches were started at 19.4 MeV. In this part of 


the analysis, due to the large amount of available data, only 
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42 
energies for which both differential cross section and polarization 


data were available were included in the analysis. This is important 

since the best differential cross section and polarization fits do 

not always have compatible parameters. The optical model parameters 

resulting from these analyses are given in Table 4. The corresponding 

best fit curves are shown in Figures 9 through 17. The curves for 

pt She at 19.4 MeV, 30 MeV, and 156 MeV will be presented in Chapter 6. 
The curves on the whole fit the experimental data quite well, 

The differential cross section fits are generally too small in the 

backward angles. The polarization fits are excellent with the single 


exception of the pt She curve at 49 MeV. It is noted that the 49 MeV 


pt She data was extremely difficult to fit using reasonable parameters. 
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TABLE 4 43 
OPTICAL MODEL PARAMETERS 


Rar snes tm ike 1.501 fm e 1.049 fm roar hes: fin 


fe) SO 
a, = 0.144 fim a. = 0.378 fm a, = 0.289 fim 
EY Yo My Wp VS0 5) % 
(MeV) (MeV) (MeV). = (MeV).=—Ss (May) No Np 
for +h 
6.0 50255 ORo2 0.0 [aes ay = 
o...0 47. 36 1.49 0.0 2.47 6.8 8 
18.0 42.98 4.69 0.0 ey (SE = 
TOW5 42.07 5.46 0.0 3.07 | - 
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Figure 9 Optical model fit to polarization data for n +34 at 
22.1 MeV. Data (Se72). 
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Figure 10 Optical model fits to differential] cross 
Data (Se72). 
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Figure 11 Optical model fit to 5.51 MeV p He scattering 
data (Mc64, C164, Mo69). 
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Figure 12 Optical modél fit to 6.82 MeV p +He scattering 
data (Mc64, M069). 
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Figure 13 Optical model fit to 8.82 MeV p +"He scattering 
data (Mc64,Mo69) . 


ay gee tt ean mn renee 
A Hi) ; 1 Nee OMe t } 


Sa ae ee ethane emp samtee tienes wee ar iste So tteeetors rar 


49 


lO.77 MeV. 
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Figure 14 Optical model fit to 10.77 MeV p +He scattering 
data (Mc64, Mo69). 
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Figure 15 Optical model fit to 13.6 MeV p'+°He scattering data 
(Hu71, Ti68). 
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Figure 16 Optical model fit to 16.2 MeV p +e scattering data 
(Hu71, Ti68). 


Dc 


100. 


(mb/sr) 


POLARIZATION 


O Oa I20 ISO 
CM 


Figure 17 Optical model fit to 49 MeV p +°He scattering 
data (Ha70, Mo70). 
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CHAPTER 6 
EXCHANGE 


The standard optical model cannot easily explain the large 
backward scattering cross sections observed. This is clear from the 
figures presented in Chapter 5. From the resonating group studies 
we find that when the wave function is antisymmetrized; the large 
back angle scattering can be reproduced. Thompson and Tang (Th71) 
found that the resonating group formalism leads to a differential 
equation containing a direct scattering potential as well as three 
additional terms. These three terms correspond to (a) nucleon exchange 
(knock out) (b) heavy particle pick up and (c) "nucleon-rearrangement" 
(see Figure 18). Using the first Born approximation they found that 
(a) is forward peaked and (b) and (c) are backward peaked. These can 
be replaced in the Born approximation by equivalent energy dependent 
potentials. Thus, (b) and (c) give rise to a Majorana exchange with 
a Genny dependence. 

Greenlees and Tang (Gr71) performed the following calculation 
to obtain a form suitable to use in an optical model analysis. Recall 
from resonating group studies: 


2 
a Ve +E- Vi(r) - Ve(r)}y(r) = [k(ror' )y(r)dr' (6.1) 


where yp is the reduced mass of the particles, E is the relative 
energy of the nuclei in the center of mass frame; ve is the coulomb 
potential; Vp is the direct potential, and K(r,r') is the kernel 


function which contains additive terms of the type: 
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n+@ SCATTERING 


DIRECT PROCESS 
(V51, "52, %53, “5q) 


KNOCKOUT PROCESS 
(v5) 


HEAVY —- PARTICLE PICKUP PROCESS 


‘52, "53, "54 ) 


“i2, “is, “ie 
Figure 18 Schematic representation of then +a 
scattering problem as a sum of direct 
and exchange processes. The nucleon- 


rearrangement process is not shown. 
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3 9) (rev! exp [-B, (r-r!)*] (6.2) 
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go(rsr' Jexpl-B,(r+r!)*] (6.3) 


ras 
ine) 
— 
s 
. 
> 
— 
il 


where 94 and Jo are form factors symmetric in r and r', and By and 
Bo are non-local range parameters. 

Since the antisymmetrization procedures may be important for 
heavier systems it is reasonable to assume that the general non-local 
optical potential should contain a direct potential Vp as well as a 


kernel of the form: 


Keb pak (ror) + iar) (6.4) 
oS -yr!/)2 -B O72 
= g,(rar')e id a : Av avlasin ASS (6.5) 


where the second term has not been previously considered in non-local 
optical models, it is particularly important in scattering at backward 
angles. 

Let fos fy and fy be the scattering amplitudes corresponding 
to the potentials Vp» Ky» and K,. At high enough energies, one may 


use the Born approximation to calculate these: 


Th aeue Ey, [exp (-ikg-r)Vp(rexp ik, “rar (6.6) 
em 
f, = - —§ lexp(-ik,-r)K,(r.r' )exp(ik.-r)dr dr' (6.7) 
] 2 he al ee mil eater 
27m 
fet exp(ike:r)K,(r,r')exp(ik.-r)dr dr' (6.8) 
2 ont Ab re lie oe Sy et, eee 


where kK.» and Ke are the initial and final wave vectors. 
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Now fy is large only in the forward direction if Vp is smoothly 


varying, and f, and f, can be written: 


oer tog Jexpty (key) (et-0) 


exply (Ky-ke)(r'tr)IKy p(rer' Jar drt (6.9) 


At forward angles where k. Ss K the contributions to fy and fo come 
mainly from the integration region where r = r'; therefore fy is 
important at forward angles and f, is not, since Ky >> K,. At 
backward angles where k. x -ke the contributions to the scattering 
amplitudes are mainly from the integration region r ~ -r' where 

K 


<< Ko. Hence, f, is important at backward angles and fy 15, tnot. 


] 
As an example, choose: 


Vp = “Vg exp [-op 7] (6.10) 
J, = 949 exPl-a, (r@#r'*)1 (6.11) 
3 (hy es 
Go = -Goq exP[-co(rotr'”)] (6.12) 
these yield: 
3/2 peek)! 

fy baiae V0 az) exp eee (6713) 
n2 3/2 -(ketky) kek)? 

f, = 9, ,( ) exp| | exp | 

] Oath 10 ay (a, +28 ) Ba, + 168, 0, 


oan a a ane a ne 
a¢ ne me me a | 
s fonts at Y tt ort birt ats tpsom ow ra 


ae e he PS t eg 


(2.9) 


eish n ah ean ea 
eh yt mor analtt oY apt 0 heate 
th. << OO iRane lahat 26 a 
ge teetese oid! of din dated 4 Freate 2s 
evar ‘n+ = 9 GOs wGhden ei | 


dona? , > bad aehpele Rene moh i om E 
Fei on wd 


(€1.@) e.g ‘se ita 2 By 


2 B/2 al kgtk, )? -(kerk, )? 
f, = donk a )  exp|—— exp |g———+ 
2 = 542 920 ‘ata, 7B) “say | ©? Laas = Ten, 
(6.15) 


verifying the above conclusion, (Note that the resonating group 
studies indicate that 91 “Qo and By = Bo). 

Therefore, it is important to include a Ky type kernel in | 
the non-local optical potential to account for large backward angle 
cross sections. 

For local potentials, we note, in the Born approximation, 
that f, and f, in Equations (6.14) and (6.15) are also the scatter- 


ing amplitudes of the potentials: 


to? alieck® Oy 2 
let) it “0G exp|- 2 on *2By z|er(- By (0, +28, )r | (6.16) 
3/2 2 0, 
i T Taertk pte VANE 
Vo(r) = “92 (30) exp|- at) exp| ay 2282) Jp (6.17) 
where k = |k.| =|ke| and P’ is the Majorana operator (space exchange). 


This suggests that, it may be possible to simulate the exchange 
effect by allowing the optical potential to contain a Majorana 


component. Thus we may write: 


V=Vi(r) + V(r) pene (6.18) 


A simpler, but cruder alternative is: 
V = (a +bPY)Vp(r) (6.19) 


where Me and Vi are smoothly varying potentials, with a and b as 


adjustable energy dependent parameters. 
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Clearly, the addition of a term of (ye dependence can 


produce a rise in the back angle cross section since it can reduce 
the cancellation of the alternate Legendre polynomials in the 
backward angles. 

In the present work, the exchange terms were included in 


the potentials in the following way: 


QL 
Mie “+ Vo e + Cail) ] 


QL 
es WOT + C{-1) ] (6.20) 


Q 
Y50 ie Veo ll ¥ Coote ] 


The exchange term is included in the imaginary component since it 
has been found necessary, by other authors (Br72, Th73) to obtain 
a good fit. These terms were incorporated into the computer code 
TASHTASH. 
The n+ 3H case was not analysed with the exchange terms 
Since the data does not extend into the backward angles, sufficiently. 


She data sets. Sig- 


A manual search was performed over all the p+ 
nificant improvements were observed for only two energies, 30 MeV 
and 156 MeV, which are shown in Figures 19 and 20. Improvements for 
the other energies were not as impressive. The case of the 19.4 MeV 
curve is shown in Figure 21. The values obtained for the exchange 
parameters are given in Table 5. Figures 22 through 27 show the 
individual effects of C, Cy and Coop for 30 MeV and 156 MeV in 
comparison with the pure optical model fits. (Note that the optical 


model parameters did not change when searched upon once the exchange 


parameters had been determined. ) 
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TABLE 5 
EXCHANGE PARAMETERS FOR P +2HE 


ED CG Cy CC ey Sp dy 
(MeV ) (MeV ) (MeV) (MeV ) Oo p 
19.4 0.0 0.0 O70 SUS. 248.8 
19.4 -0.012 -0.25 0.0 340.8 61.3 
30.0 -0.065 -0.4 0.06 25.2 18 


156 -0.066 0.007 0.023 74.5 24.9 
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FIGURE 19: Optical model fits for p He at 30 MeV with and without 


exchange. Data (Ha70, Mo70, Ha68). 
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Figure 20 Optical model fits for p +3He at 156 MeV with and 


without exchange. Data (La70, Frpc). 
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Figure 21 Optical model fits for p + He at 19.4 MeV with and without 


exchange. Data (Ba71). 
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FIGURE 22: Optical model fits for p +*He at 30 MeV with 


Ce -0.065 MeV and without exchange. 
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FIGURE 23: Optical model fits for p He at 30 MeV with 


Cy = -0.4 MeV and without exchanae. 
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FIGURE 24: Optical model fits for p He at 30 MeV with 


ee = 0.06 MeV and without exchange. 
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Optical model fit for p +*He at 156 MeV with 


Cy 


= 0.004 MeV, 
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FIGURE 27: Optical model fit for p +"He at 156 MeV with 


ba = 0.018 MeV. 
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In the 30 MeV and 156 MeV cases, a marked rise in the 
differential cross section at backward angles is observed when 
the exchange is included. For the 30 MeV polarization, the 
exchange improved the fit tremendously. For the 19.4 MeV, the 
exchange improved the fit to the polarization, however, the fit 
to the differential cross section was made worse. For energies 
where this disagreement occured, the exchange yielded no overall 
improvement and was, therefore, left out in the best fit. From 
Figures 22 through 27, one sees that the Co and Cy were the major 


contributors in the backward angle cross section increases. Also 
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these were the major contributors in improving the 30 MeV polarization 


fit. 
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CHAPTER 7 
DISCUSSION AND CONCLUSION 


DISCUSSION 

In this section, comparisons will be made between the curves 
obtained by the various analysis methods and the optical model. 
Exchange in the optical model will also be considered. 


n +34 


All the discussed analysis methods were performed for this 
system on the differential cross section at 6 MeV. The optical model 
curve is one of the better fits to the data. All the available 
differential cross section data for this system have been analysed by 
LeMere et al (Le75) and ourselves (Sh72). Over the whole data set the 
optical model fits are at least as good as the resonating group fits. 

We also analysed the only polarization data available, at 22.1 MeV, and 
achieved a good fit. 

One must note that no data sets were available for which 
there was polarization as well as differential cross section data. This 
does allow better fits since in many cases the polarization data requires 
the parameters to go in a different direction than the cross section data 
requires. This is not too serious a problem in this case since the 
parameters obtained at 21 MeV are not very much different from those 
obtained at 21 MeV and 23 MeV (see Tabte 4). 

On the whole, the optical model gave very satisfactory fits to 


the available experimental data from 6 MeV to 23 MeV, for n +34 elastic 


scattering. 
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For the p +3He system, there exists a great deal of data, 
therefore, the analysis was restricted to energies for which both 
polarization and differential cross section data were available, i.e. 
from 5.51 MeV to 156 MeV. Note that the well geometry used is the same 
as for the n 4 system. 

The agreement between the model predictions and the data is 
good, with the exception of 49 MeV, for which a good fit could not be 
obtained with reasonable parameters. Inspection of Table 4 (for both 
p +3He and n +3H) reveals the following trends in the potential depths: 

(1) The real well depth decreased uniformly as the energy 
increases. 

(2) The imaginary volume well depth, Wy increased in general, 
as expected. 

(3) The imaginary surface well depth, Wy» Was kept at zero for 
the best fit. 

(4) The spin-orbit well depth increased smoothly with energy. 
A quadratic least squares fit was performed on the variation 


of the real well depth with energy. The resulting equation is: 
V,(E) = 53.75 - 0.869 E + 0.004 £2 (7.1) 


This equation can be compared with those obtained by other authors, for 
proton scattering from: 
tHe (Sa68): V(E) = 53.3 - 0.69 E (7.2) 
Medium-heavy nuclei (Ro66): ve) = 53.8 - 0.33 E (7.3) 
It appears that as the mass of the target nucleus decreases, the slope 


of the curve becomes steeper. Also, note that if a linear fit was used 
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72 
for the p + He case, the slope would have been somewhat flatter (and the 


fit not nearly as good). See Figure 28 for the plot of Equation (741 )e 

The exchange effect made a marked improvement in the 30 MeV 
and 156 MeV optical model fits. The fits obtained including non-zero 
exchange terms for the other energies did not really improve the curves. 

The 49 MeV case has presented many problems since it refuses 
to be fit decenty. This leads one to believe that the data maybe 
incorrect, however the 48.5 MeV differential cross section data agrees 
almost exactly with the 49.5 MeV data. Note that in the 49 MeV analysis, 
49.0 MeV polarization data was analysed with 49.5 MeV differential cross 
section data. It is not at all clear what is the problem at this 
energy. 

Comparing the microscopic calculation fits for polarization 
and differential cross section for 5.51 MeV to 10.77 MeV with the optical 
model fits, one finds that the microscopic calculation fits are better 
at lower energies and equivalent for the higher energies. The optical 
model does have an advantage since it is easily extended to higher 
energies where the microscopic calculations are not. The optical 
model fits are superior to the resonating group fits for energies less 
than 12 MeV, with the advantage that it can handle the polarization data 
that the resonating group cannot. Comparing the phase shift analysis of 
p + He at 30 MeV with the optical model analysis, one finds that the 
phase shift analysis obtained five equivalent type 1 fits, all of which 
were better than the optical model curves (with and without exchange). 
However, the phase shift method required nineteen parameters and obtained 
at least four extraneous solutions which cannot be eliminated. Whereas, 


the optical model fits are good and the problem of extraneous solutions 
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Figure 28 Least squares fit to the V5 8 obtained from the 
model fit for p +He; yielding the equation: 
V(E) = 53.75 = 0.869 E + 0.004 E° 
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| does not exist. 


CONCLUSION 

This work was undertaken in order to demonstrate the usefulness 
of the optical model as an analysis method for scattering off very light 
nuclei. We have found that it yields fits which are comparable or 
Superior to the standard methods. The optical model is more useful than 
the resonating group method (also microscopic calculations) in analysis 
over a broad energy range. It is easier to use than phase shift 
analysis since it does not employ as many parameters or obtain 
extraneous solutions. 

We believe that the optical model is a valid semi-phenomenological 
analysis method for the elastic scattering of nucleons off very light 
nuclei, and when coupled with exchange terms becomes a very useful tool 
indeed. 

A logical extention to this work would be to include spin-spin 


terms in the potential in order to determine their importance. 
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APPENDIX 
ALTERNATE DERIVATION OF DIFFERENTIAL CROSS SECTION 


The following derivation of the differential cross section 
is based on the evaporation model, as described by Vogt (Vo68). 
Define a wave function bo» corresponding to a channel c: 

We = Me oF YT Cesmm jum) yx, (A.1) 

uy mM. ie Ss 
where re is the relative separation of the pair of particles in the 
channel c; vy is the state of internal motion of the two particles 
in the channel; Yn, is a spherical harmonic; Xsme is the wave 
function of the coupled intrinsic spins 1 and i of the two parti- 
cles, and (gsm,m, | dm) is a Clebsch-Gordon coefficient. 

Also, the total wave function of the system can be written: 


Us Ih Cee (A.2) 


Cc 


where Y. are the radial wave functions for each channel. 
Considering only coulomb and angular momentum barriers, the radial 


wave equation for channel c is: 


2 2n..k 
Pe. r ur Cis 
Cc Cc 
where 
Zig ef 
odes ras 
Ne Fv ¢ 


is the usual coulomb parameter. 
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80 
Equation (A,3) has the regular solution Fe and the irregular solution 


Go. These solutions can be combined to give the incoming wave I. and 


the outgoing wave 00: 


iw. 
ls =!) b= (G. = iE So) e (A.4) 
where 
Q (n=) 
Cc - 
ve = tan 1s (A.5) 


and is the coulomb phase shift. The incoing and outgoing waves have 


the asymptotic form: 


pS wal: Z 
Lies 0. = exp I sie = On In 2k or) ; (A.6) 


Since each ie can be expressed as a linear combination of 


incoming and outgoing waves, we can write: 


] 
pee MY once tava hc (A.7) 


where A. and B. are amplitude coefficients and the factor ve 


insures 
unit incoming or outgoing flux. Solution of the Schrodinger equation 


for the whole system yield all the B.S» Tec: 


= ) Un AY (A.8) 
c! 


where Une: is the collision matrix. 
Consider a nuclear reaction which has an initial pair of 
particles a and a final pair of particles a' observed in a direction 


(0, ¥) relative to the incident beam. The total wave function can 


be written: 


. PaT 
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where ee is the incident plane wave and 
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In order to obtain the Aa Ss Equation (A.10) is compared with 
the incoming and outgoing waves at distances where the coulomb fields 


of the target and projectile are screened from each other. This gives 
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Substituting Equation (A.11) into Equation (A.10), we obtain: 
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for the a, s, and m. of the incident beam; A. = 0 for all other channels. 


Also, 


J+s ie 


J 
: -] 35 
Bow = Busepy talked” (avg)? T (280m, [dmy) (2241 PU ey. gage - 


(A.14) 


ig ; 4 


(@a)- ‘i | 


” ‘ sit J On) ad See 


(Of A} 


atte besgqmes ef (or. A) not aup8 o2 A. os aera “ "1 
ebtat? dnofvos sit svat esondtel> ‘6 davis antopsuD ? 
v930 rope movt benesio2 sis sthiesterg:t 


eevip erdl 


~ 


(Tf.A) 
-atagda om (OFA) nova oon en 


5" as tha) ed Yale 


(SEA) 


7 aA r fo 


“Chas signi goo, wie * we ae 7 tie | 


Pita 


iB A) 


.éfannata yertte Tis vet 0 = of nod taebtont al? to 0 elon he yi. orit yo? 
; ; ; a 


Tae 
Pee 


ee) eins oe iy v 
» 
4 ' ir % f] oa 
1 ) y ‘4 . 
"es j i i | ay ey mer as a 
i ae! a eS 


q 
is 
ub) 
a” 
yee 
‘4 te ' 
i 1 Bi 
vey 
th 
pS 
“4 J 
' 
es, | 
‘ 
t 
a at 
‘ LL 


82 
Using Equations (A.9), (A.12) and (A.14), we can write: 


(6 aU eeween ) (A.15) 
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The kronecker delta accounts for outgoing waves in the incident plane 


wave. A more useful form for w is in terms of the laboratory 


reac 
variables: 
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Now, the differential cross section for the reaction o > a' 


can be written: 
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This equation can be simplified by performing some of the sums over 


the Clebsch-Gordon coefficients, to obtain: 
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where Z's are the Z coefficients defined by Biedenharn, Blatt and 


Rose (Bi52) 
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For elastic scattering: 
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Tombrello obtained an equation for the differential cross 
section for spin-1/2-spin-1/2 scattering in the center of mass 


system (To62): 
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where Se s'32,s) is the element of the scattering matrix for total 


angular momentum j connecting initial orbital angular momentum &' 


and channel spin s', with their final values 2 and s. In these 


formulae 
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where v is the relative velocity of the two particles, and wu is the 
reduced mass of the system . 


The polarization of the scattered protons (To65) is given 


by 
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